Gross features of finite nuclei at finite temperatures 
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Abstract 

A simple expression is obtained for the low temperature behavior of the energy and entropy of 
finite nuclei for 20 < A < 250. The dependence on A of these quantities is for the most part due 
to the presence of the asymmetry energy. 
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In nuclei the nuclear force between the nucleons is short-ranged which leads to saturation 
of the binding energy, E{A,Z), per nucleon. Empirically it is known that for the stable 
isotopes 

(1) 

where a ~ 8 MeV and A > 20. This fundamental gross property of nuclei is well accounted 
for by the semi-empirical mass formula [H El [3] which provides a simple parametrization of 
the binding energy per nucleon for all known nuclei. Qualitatively it is also consistent with 
the simple Fermi gas model prediction [U [3] 

where S is the ground state energy and e/ is the Fermi energy which is constant as long 
as the particle density remains constant. It is interesting to note that at low temperatures 
{T < Tf) the excitation energy in the Fermi gas model is given by [U E] 

S{A,T) 3 TT^T^ 
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or (3) 



= a + bT^ (4) 

where a and b are constant again as long as the particle density remains constant. Although 
the Fermi gas model may be an oversimplified model, nonetheless it underscores the relevance 
of independent particle (or quasi-particle) methods in nuclear structure physics. Mean field 
methods have been used throughout the periodic table, both at zero temperature and at 
finite but low temperatures, and typically yield a behaviour of the energy density. This 
suggests that perhaps at low but finite temperatures a simple scaling relation might also 
exist for the excitation energy, S{A,T), of finite nuclei. 

In order to test the validity of equation Q we have made use of a finite temperature 
extension of the semi-empirical mass formula [7]. The following form for the temperature 
dependent binding energy [3l E] has been assumed 
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EiA,Z,T) = a(T)A + /5(T)A3 + (7(T)-^)( 



r^{T)^Atl+A\t^\^ 
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+ K(T)-,{1 - - • ; ) + 5(T)/(A, Z)A-^ (5) 

As Zi (0.8076)2^3 
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where A = N + Z , = ^{Z — N) and f{A, Z) = (—1, 0, +1) for (even-even, even-odd, odd- 
odd) nuclei. Here 1 is the volume energy, 2 is the surface energy, 3 is the asymmetry energy, 
4 is the Coulomb energy and 5 is the pairing energy contribution to temperature dependent 
binding energy. In this parametrization the temperature dependence of the contributions to 
the Coulomb energy term which arise from exchange and surface effects [31 19] was ignored. 
Also, no attempt has been made to include shell effects in the finite temperature expression. 
The excitation energy per particle is given by 

£iA Z,T) ^ EiA,Z,T)-EiA,Z,0) 

A A ' ^ ^ 

At T = the coefficients are given by [3J a(0) = -16.11 MeV, (3{0) = 20.21 MeV, 
7(0) = 20.65 MeV, r]{0) = 48.00 MeV, and k{0) = 0.8076 MeV obtained from a fit to the 
experimental nuclear ground state energies of 488 odd mass nuclei. The T = coefficient 
for the pairing term is taken as 5(0) = 33.0 MeV [8]. 

To obtain the temperature dependence of the coefficients, the available experimental 
information about the excited states of nuclei throughout the periodic table was used to 
determine the partition function of each nucleus in the canonical ensemble 

Z{A,Z,T) = Y,9^^M-^^E,) + dE gA,z{E)exp{-(3E) (7) 

i 

where Qi = 2ji + 1 is the spin degeneracy factor and Ei the excitation energy of the ith 
state of the nucleus, and /? = 1/T. The sum in the first term of equation ([T]) runs over the 
experimentally measured (discrete) excited states. 

Since the experimentally known spectrum in most cases is only sufficient to allow the 
accurate determination of Z for very low temperatures (T <^ 1 MeV), it is necessary to 
supplement the experimentally known spectrum with an appropriate approximation to the 
continuum g^ ziE). For this purpose, the fits obtained in [TU] were used. For sufficiently 
large energies, the usual Fermi gas expression for the total density of states (i.e. including 
the spin degeneracy) is used: 

v/^ exp(2^a^,^[/) 

Here aA,z is the level density parameter and U = E — P{N) — P{Z), where P{N) and P{Z) 
are the pairing corrections for neutron number and proton number Z respectively [TO] . 
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This parametrization is obtained by means of a saddle point approximation, and is probably 
only valid up to T ^ 6 MeV. However, in the region up to T ^ 4 MeV, this parametrization 
is probably acceptable. 

At lower energies, a suitable fit to the nuclear energy level density can be obtained to the 
form 

v^g exp {E - Eq) 
9A,z{E) = , (9) 



T T 

with a the spin-dependence parameter. Values for the parameters aA,z, Eq and a, as well 
as the respective regions where ([s]) and ^ should be used, for a large number of nuclei can 
be found in |10] . 

The nuclei used to determine the aforementioned coefficients can therefore be divided into 
three groups: Nuclei where sufficient discrete states are known to allow the use of the discrete 
spectrum at low energies and the Fermi gas expression ^ at higher energies. Nuclei where 
the discrete spectrum does not extend high enough for ^ to be valid. For these nuclei, the 
discrete spectrum is used for low excitation energies, followed by the exponential form ^ 
for intermediate energies and finally the Fermi gas expression ^ at high energies. Nuclei 
where very little of the discrete spectrum is known. In these cases, (joj) is used for the low- 
and intermediate- excitation portions of the spectra and ([s]) for the highly excited part. All 
three groups are spread across the whole periodic table. 

The lower bound En on the integral 
in ([T]) is taken to be the energy at which 
(|8| should become valid (from [lOj) for 
the first case above, 80% of the largest 
discrete energy level for the second case 
above, and zero for the third case. For 
temperatures up to T ^ 4 MeV, the 
upper bound E^ax ~ 3 GeV was used. 
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The coefficients in the mass formula FIG. 1: Quadratic fit of the calculated binding energy of 
have been determined by a least squares xe (A=130) using equation (§. The fit to the Xe data is 
fit of (§ to the ensemble average of the representative of the fits obtained for other nuclei. 
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determined from a total of 313 nuclei in the mass region 22 < A < 250 for temperatures 
T < 4 MeV. The temperature dependence of the six coefficients is given in reference [7j. This 
finite temperature parametrization has been used to identify the pairing phase transition in 
symmetric nuclear matter [11]. The same techniques have also been successfully employed 
to identify the remnants of the pairing phase transition in finite nuclei [12] • 

The excitation energy, S{A, Z,T), has been calculated for a number of stable isotopes 
in the mass region 20 < A < 250 using the temperature dependent coefficients determined 
in reference [7]. For each stable isotope, the constant b in equation ^ was determined 
from a quadratic fit of S{A, Z,T) versus T. These fits were in general excellent, as shown 
in figure [T] with only occasional small deviations from the quadratic fit occurring at low 
temperature. This is not surprising given that the continuum contributions are given by 
Fermi gas expressions. However, this alone in no way guarantees any simple dependence of b 
on A. Note also that no attempt has been made in the fits to take into account the presence 
of low temperature collective to non-collective phase transitions fn\ [T^ . The coefficient b 
was then plotted as a function of A, as shown in figure [2j As can be seen in figures |2] and [3] 
there is a relatively simple dependence of b on A. In order to determine the effect on b of 
the different terms in the temperature dependent binding energy ([s]), various combinations 
of the five terms were plotted against T then quadratically fitted to obtain b for each A and 
finally these points were linearly fitted, as shown in figures [2] and [3j Although only 20 stable 
nuclei were used these were chosen at random to cover the range from ^4=20 to A=250. 
In figure [3] the equations for the fits are b = —2.032 x lO^^A + 0.152 when all terms of 
equation ([s]) are used and b = —1.617 x 10~^A + 0.140 if term 3 of equation ([s]) is excluded. 
The gradient is small in both cases but especially small when term 3 is excluded, with b 
gradually varying from 0.141 to 0.137 for nuclei larger than ^4=40. Surface effects become 
significant in nuclei smaller than A=40 and the pairing effects are most likely overestimated. 
This most probably accounts for the deviation of the corresponding points from the fit in 
figure [3] and subsequently prompted the exclusion of these points from the linear fit. These 
results demonstrate that the asymmetry term (i.e. term 3) is largely responsible for the A 
dependence of b. 
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FIG. 2: The resultant b from various combina- 
tions of terms in equation (|5| has been plotted 
against A, thereby illustrating the linearity of all 
five cases. 
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FIG. 3: Linear fit of the b resulting from 
equation (|5| and equation (|5| without term 3 has 
been plotted. 



Furthermore in the Canonical ensemble the specific heat is given by 

c - I (U) 

At low T up to an additive constant the entropy per particle, S, for finite nuclei in the mass 
range 20 < A < 250 is therefore given simply by 

S = \hT (13) 

with h = -2.032 x 10"^ + 0.152. 

It is interesting to note that the simple Fermi gas model does not take properly into 
account that the nucleus is composed of protons and neutrons and that therefore the effects 
of the Coulomb force must be considered. Were only the nuclear force present one might 
expect b to be independent of A. The presence of the Coulomb force gives rise to the 
asymmetry term which at T=0 leads to a displacement of the island of stability away from 
nuclei with N = Z and at finite temperature is largely responsible for an A dependence in 
b. 
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